Nonlinear amplification of instabilities with longitudinal expansion 
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We study the dynamics of nonequilibrium instabilities in anisotropically expanding systems. The 
most prominent example of such a system is the 'Glasma' in the context of relativistic heavy-ion 
collision experiments, where the expansion is a consequence of approximately boost-invariant initial 
conditions. Here we consider the problem of parametric resonance in scalar A'^-component quantum 
field theories with boost-invariant initial conditions, which is similar in spirit. We find that many 
aspects of the dynamics can be treated analytically by introducing a generalized conformal time. 
Primary instabilities, which are described by the linearized evolution equations, are seen to lead to a 
secondary regime of amplifications with strongly enhanced growth rates due to nonlinear corrections. 
For the secondary instabilities we present a power-counting scheme for weak coupling, and discuss 
their role for the question of isotropization and the establishment of an equation of state. 



I. INTRODUCTION 

Plasma instabilities in quantum chromodynamics 
(QCD) can play an important role for our understanding 
of the early time, and maybe even rather late time, evolu- 
tion of matter created in ultrarelativistic heavy-ion colli- 
sions at sufficiently high energies [1]. Quantitative stud- 
ies of the evolution at early times are typically performed 
in either numerical simulations of the classical-statistical 
field theory [2H5] or in simulations of semi-classical trans- 
port approaches such as the hard-loop framework [6J and 
taking into account backreactions on the momentum dis- 
tribution of the hard particles [7] . 

A crucial open question in this context concerns the 
role of quantum fluctuations during the early far-froni- 
equilibrium evolution of the collision. While the origin of 
an instability can typically be understood from unstable 
modes of the linearized evolution equations, the inclusion 
of nonlinear quantum corrections can lead to important 
phenomena. Most strikingly, the primary growth of fluc- 
tuations from unstable modes can lead to a secondary 
stage of strongly enhanced amplification in a wide mo- 
mentum range. This can have important consequences 
for the question of isotropization and the establishment 
of an equation of state in the context of heavy-ion colli- 
sions. 

The phenomenon of secondary amplification from non- 
linear quantum corrections has been studied in great de- 
tail for scalar quantum field theories in the context of 
early- universe infiaton preheating [8j |9]. In that case 
the dynamics with isotropic expansion can be mapped 
onto an equivalent problem in Minkowski space-time with 
the introduction of a conformal time for massless the- 
ories. Nonlinear amplifications with strongly enhanced 
growth rates have also been observed in fixed-box stud- 



ies of plasma instabilities in SU{2) and SU{3) Yang- 
Mills theory [31 [TU] , where secondary growth-rates have 
been pointed out as a possible mechanism to speed-up 
isotropization. Similar phenomena may also be identi- 
fied taking into account the longitudinal expansion asso- 
ciated with approximately boost-invariant initial condi- 
tions in classical simulations [H |5] . Here the longitudinal 
expansion of the plasma leads to a dilution of the system 
which naturally competes with the instability. Hence it 
is non-trivial whether primary instabilities can produce 
high enough occupation numbers for secondary processes 
to be relevant. 

In this work we investigate generic phenomena that 
occur in anisotropically expanding systems undergoing 
an instability. We consider scalar A^-component field 
theories with parametric resonance initial conditions, for 
which many aspects can be understood analytically from 
the underlying quantum field theory by introducing a 
generalized conformal time. Primary instabilities, which 
are described by the linearized evolution equations, are 
described in Sec. |Tl] An important new ingredient go- 
ing beyond this linear regime is the appearance of the 
time-dependent spectral function as an additional lin- 
early independent correlation function in the evolution 
equations. In Sec. [TIT] we present an analytic discussion of 
the nonlinear amplification of instabilities, where the pri- 
mary growth of fluctuations is seen to lead to secondary 
instabilities with strongly enhanced growth rates. We 
present a power-counting for weak coupling, which takes 
into account the time dependence of the nonlinear cor- 
rections. 

The nonlinear amplification is followed by a nonper- 
turbative regime, where no power-counting for weak cou- 
pling can be given. For the scalar theory this regime can 
be described within a 1 /N expansion of the two-particle 



2 



irreducible (2PI) effective action [TT] to next-to-leading 
order }12) . The dynamics beyond the linear regime can 
also be described using classical-statistical simulations, 
which are also available for gauge theories [2H5]. For 
comparison with the analytic estimates, we employ a 
classical-statistical approach for our numerical simula- 
tions in Sec. lIVI We conclude with Sec.|Vl Further details 
about the generalized conformal time in anisotropically 
expanding coordinates as well as the analytic solutions 
are given in appendices A to D. 



II. PARAMETRIC RESONANCE WITH 
ANISOTROPIC EXPANSION 



A. Scalar field theory and initial conditions 

We consider a real-valued, A^-component scalar quan- 
tum field theory with quartic self-interaction. The clas- 
sical action is given by 



(1) 

with self-coupling A for the fields Lpa{x) depending on 
space-time coordinates = {t, x) and components a — 
l,...,iV. Here g^i,{x) is the metric tensor and its de- 
terminant is denoted by g{x) = dct gfj,^{x). Since boost- 
invariant initial conditions along the longitudinal direc- 
tion z = x^ will be of interest, it is convenient to intro- 
duce the comoving coordinates 



rj — artanh ^-^ 



(2) 



The fields becomes then a function of 'proper time' 
T, rapidity rj and the remaining transverse coordinates 
XT — {x^,x'^). The metric tensor in comoving coordi- 
nates follows from the transformation from Minkowski 
space as g^^^ = diag(l, —1, —1, — r^) and ^/— g = r. 

We are interested in the real-time evolution of the cor- 
responding quantum field theory for scalar A^-component 
Heisenberg field operators ^^(t, xt, '7)- For a given den- 
sity matrix pd{tq) at initial proper time tq the field ex- 
pectation value is {(I){t, xttIj)) = Tr[p£)(ro)0(T, x^, »7)] 
and equivalently for products of field operators, which 
determine the correlation functions. We will consider 
To > 0, i.e. the time-evolution in the forward light-cone, 
and spatially homogeneous expectation values. By spa- 
tially homogeneous we mean homogeneous in the trans- 
verse plane and homogeneous in longitudinal rapidity at 
given proper time. We employ Gaussian initial condi- 
tions, which can be conveniently formulated in terms of 
the macroscopic field 



(3) 



as well as the statistical two-point correlation function 
Fabir, t', XT - x'rp,-!]- r]') 

= I ({Mr,XT,v),MT',x'T,r^')}) - Mr)Mr') (4) 

and derivatives at initial time. Here {.,.} denotes the 
anti-commutator. Using 0{N) symmetry in field-index 
space, we take the macroscopic field to point in the a = 1 
direction, i.e. 



(5) 



The fluctuations can be taken to be diagonal, i.e. Fab — 
diag(Fj|, . . . , F±) with the subscripts, longitudinal 
and transverse, indicating the orientation with respect 
to the macroscopic field. 

Relevant for parametric resonance are weak couplings 
(A ^ 1) with a parametrically large macroscopic field 
(02(ro) ~ C(l/A)) and smah fluctuations (F(ro) - 
at initial time. For the field initial conditions we write 



dr(, 




'/'(t'o) 



(6) 



where the derivative at r = tq mimics an initially free- 
streaming behavior and uo denotes the rescaled initial 
field amplitude. If not stated otherwise we will employ 
the typical set of parameters = 0, (JoTq = 5, N = i 
and A = 10^^ for our estimates and numerical simula- 
tions below. The spectral shape of the initial correla- 
tion functions will be of little relevance as long as the 
magnitude is sufficiently small compared to the initial 
macroscopic field squared. For completeness the spec- 
tral distribution of the considered fluctuations is given in 
appendix D. 



B. Linear regime and generalized conformal time 

Since the initial macroscopic field is large and fluctu- 
ations are small, at sufficiently early times the time evo- 
lution is accurately described by the classical evolution of 
the field and linearized fluctuations around it. The sub- 
sequent onset of non-linearities and impact of quantum 
fluctuations will be discussed further below. In the lin- 
ear regime the evolution of the macroscopic field follows 
directly from the stationarity of the classical action ([T]) 
and is given in comoving coordinates by the equation of 
a damped anharmonic oscillator 
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dr 



A 

6N 



(r) = 



(7) 



The damping term characterized by a first-order r- 
derivative arises due to the longitudinal expansion and 
represents dilution of the system. In order to solve ([t]) 



3 



analytically, we perform a change in the time variable as 
well as a rescaling of the fields by introducing 



(JqcIt = 



air) 



1/3 



d9 



with the generalized scale factor 
a(T) = r 



ajr) 
a{To) 



1/3 



(8) 



(9) 



in case of a one-dimensional Bjorken expansion. This can 
be seen as a generalization of the concept of conformal 
time to anisotropically expanding systems. For the one- 
dimensional expansion the scale factor a^r) enters the 
above definition of conformal time with exponent 1/3, 
whereas in the isotropically expanding case one would 
use d/3 for d-dimensional space instead. The new di- 
mensionless time-variable 9 is explicitly given by 



2/3 



1 



(10) 



and we will refer to it as conformal time in analogy to the 
isotropic case. The evolution equation of the macroscopic 
field in terms of the new variables reads 



= 



(11) 



where we introduced the time-dependent effective 'mass 
term' 



a(0o) 



2/3 



a'{e)y ia"{e) 



a{9) 



3 a{e) 



(12) 



Here primes denote derivatives with respect to 9, which 
are explicitly given by 

a'{9) _ 1 I 29 V''^ 
a{9o) Co To \3(ToTo J 

-1/2 

(13) 



a{0o) 
a" (9) 
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a{9o) 3((7oro)2 v3c^oT"o 
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where we chose 6*0 — 3(To''o/2 in accordance with (10 1 



This change of variables is discussed in more detail in 
appendix A. 

For large initial field amplitude, (ToTq 3> 1, the above 
evolution equation can be solved appr oxim ately for the 
massless case (m^ = 0). According to (12 1 one then has 
rh^{9) = 1/(46*^), which is small compared to the macro- 
scopic field. Thus fh?{9) ~ gives a rather accurate 
description of the dynamics in this case. The evolution 
equation ( 11 ) then becomes that of an anharmonic oscil- 
lator 



di 



= 



(14) 



Its solution is given in terms of Jacobi elliptic functions, 
i.e. 




for the considered initial conditions specified in (|6| and 
appendix D. The Jacobi cosine, cn(0; a), is a doubly pe- 
riodic function in 9 with periods AK{a) and AiK{l — a), 
where K{a) is the complete elliptic integral of the first 
kindj^ Accordingly the macroscopic field ^ displays os- 
cillations with constant period in conformal time. For 
the physical interpretation of this result it is insightful 



to express the approximate solution ( 15 ) in terms of the 
original variables, where it reads 



CTO 



(7) 



1/3 



To 



2/3 



In Fig.[T]this is compared to the numerical solution of ([7| 
without further approximations. The very good agree- 
ment verifies that the terms neglected in ( 14 ) are irrele- 
vant for this choice of parameters. 

From Fig. [T] one observes that the macroscopic field 
(/)(t) decays as t~-^/^, while it displays oscillatory behav- 
ior with a constant period in conformal time 9 cx r^/'^. 
The prior follows, of course, also from simple consider- 
ations of energy-momentum conservation. The energy 
density e satisfies the equation of Bjorken hydrodynam- 



de 



Pl 



(17) 



where Pl denotes the longitudinal pressure in the local 
rest-framej^ Initially the energy density is dominated by 
the macroscopic field and proportional to 0^. As the 
macroscopic field is homogeneous and the system is con- 
formally invariant, the longitudinal pressure averages to 
1/3 of the energy density. Hence e ~ r"''/^ and, there- 
fore, 4> ^ T^^/^. The fact that the oscillations have a 
constant period in conformal time is a consequence of 
the dilution of the system. 



C. Nonequilibrium instability with expansion 

To describe the linearized evolution equations for the 
fiuctuations (|4|, we will work in Fourier space with re- 
spect to spatial coordinates, i.e. we consider the dimen- 
sionless quantity 



) = y d^XTdT^F^{T,T',XT,il) e-'(^^^^^+'^") 

(18) 

and equivalently for the longitudinal component 
f II (r, t',pt, J^)- The linearized evolution equations for 
the fluctuations can be obtained by expanding the clas- 
sical field equation of motion, given by the stationarity 



^ In terms of a series expansion, K{a) 



128 ^ 512 



(15) 

See Sec. 4 for a discussion of the stress-energy tensor. 
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FIG. 1. (color online) Solution to the linearized field equa- 
tion ^ for croTo = 5. The approximate analytical solution 
( |16[ l given by the dashed curve is practically on top of the full 
numerical one (solid curve). As a consequence of the longi- 
tudinal expansion, the field shows oscillations with constant 
period in t'^''^ while it decreases as r^^''^. 



of the action ([T]), 
homogeneous background field 



to first order in deviations from the 
r). 



Ill 



Later in Sec 

quantum corrections, these equations will also be seen 
to correspond to the evolution equations of the quantum 
field theory in the limit where all nonlinear or loop cor- 
rections are neglected. The linearized equations for the 
fluctuations read 



1 



A 



di + -dr+p'T + ^+ rn' + (r) 



5^ + Ja.+p| + ^+m^ + A^^(r) 



Fj_{t, t\pt,v) 
= 0, 

= 0. 
(19) 



The damping term characterized by a first-order r- 
derivative again appears due to the expansion. Other 
competing scales that depend differently on proper-time 
r are related to 

• the macroscopic field squared 4>^{t) ^ t^^/"^, 

• the longitudinal momentum squared jr^ and 

• the transverse 'mass' squared vr? -\-v\ which is in- 
dependent of proper-time. 

This interplay of dilution and red-shift effects induces a 
variety of new phenomena, which are not present for non- 
expanding systems. In particular, the well-known phe- 
nomenon of parametric resonance in Minkowski space- 
time will receive significant changes in the presence of 
longitudinal expansion. 

For an analytic description of the time evolution of 
fluctuations, we introduce conformal time variables simi- 



lar to what has been done for the macroscopic field above: 

where r = r(0) and accordingly r' = Ti^^')- Inserting 
the solution for the background field (15), the linearized 



evolution equations using conformal time variables read 



Fn(9.0 



\PT,v) 
= 0, 

',PT,v) 

= 0. 
(21) 



We restrict ourselves again to the massless case, m? = 0, 
and approximate fh? [0) ~ as in ( 14 ) . With this ap- 



proximation the equation of motion resembles the Jaco- 
bian form of the Lame equation. A crucial difference 
is the explicit 0-time dependence of the (dimensionless) 
momentum term 



P\0) 



29 



2 

Pt_ 

2 



4^ 



(22) 



It is instructive to consider for a moment what hap- 
pens if p^ was not depending on 0-timc. In this case (21 ) 



corresponds for vanishing •fh?[6) to Lame equations. The 
latter exhibit the well-known phenomenon of paramet- 
ric resonance, which is reviewed in appendix B. The dy- 
namics is then dominated by the transverse modes (-Fl) 
showing exponential growth as a function of time in a 
resonance band of momenta with 



< p2 



< 



1 



(23) 



For our purposes the (dimensionless) momentum- 
dependent rate for exponential growth is to good approx- 
imation given by 



l{p) 



WN) 



(24) 



where the maximum growth rate in this approximation 
is 



27re- 



70 



K{l/2) 



for 



(25) 



We give more accurate analytic expressions involving in- 
verse Jacobi elliptic and Jacobi zeta functions in ap- 
pendix B, where we also compare the approximate ex- 



pression (24) to the full growth rate in Fig. 12 



For the anisotropically expanding system the situation 
is somewhat more involved due to the effects of dilution 
and the red-shift of longitudinal momentum modes. In 
order to obtain an analytic understanding for the ex- 
panding system described by (21), it is very instructive 
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FIG. 2. Time-dependent (dimensionless) growth rate 7(pt, i^, 0) in the transverse momentum pT and rapidity wave number u 
plane for different times. The (left) panel corresponds to r/ro = 1, the (right) panel shows t/tq = 5. The instability develops 
from the high pr small v region to the high v small pr region. 



to assume that the momentum term (22 1 varies only 



slowly on the characteristic time scale of one oscillation 
of the macroscopic field. By expanding the time depen- 
dent momentum term as p^{9 -\-Te) ~ p^{0) + dep^{0)Tg, 
where Tg = AK{l/2) corresponds to one period of os- 
cillation of the macroscopic field, we find that formally 
this quasi-static approximation corresponds to the limit 
4:K{l/2)\dgp'^{8)\ P^{6). For resonant modes, where 
the characteristic momenta are parametrically given by 
Pt ^ Co B^nd V ^ cfqTq, it is straightforward to verify that 
the above condition is approximately fulfilled whenever 
cqT'o ^ 4 K{l/2), which can easily be achieved by choice 
of parameters. Within this approximation the resonance 
band and gro wth rates for the expanding system are ob- 
tained from (23) and (24) by replacing p^ -> p^{0) as 
given in (22 1. Therefore, in this approximation the time- 



dependence of the growth rate enters only through the 
explicit time-dependence of p'^{0) and is a consequence 
of the red-shift and the dilution. 



The growth rate ( 24 1 with time-dependent momentum 
terms ( 22 ) is displayed graphically in Fig. [2] as a func- 
tion of transverse momentum px and longitudinal wave 
number v for different times. One observes that the in- 
stability develops from the high px and small region 
at early times to the high and small pt region at later 
times. In the remainder of this section we will discuss 
different phenomena that appear within the quasi-static 
approach and show that these analytic results compare 
well to full numerical solutions of (21 1. 



The explicit 0-time dependence of the approximate 
resonance band criterion leads to new phenomena that 
are not present in the non-expanding case. A given 
mode with fixed transverse momentum pT and fixed 
may satisfy the condition ( 23 1 for a certain time window 
Tstart < T < TEnd, whilc thc Condition is not met out- 
side this window. Here rgtart and r^nd are related to the 



modes can shift inside the resonance band, exhibit ex- 
ponential growth for a certain time and then shift back 
out of the resonance band again so that the exponential 
growth stops. In particular, if rgtart > tq for a certain 
mode, this mode will exhibit a delay in the onset-time of 
growth. We note that this qualitative behavior is charac- 
teristic for nonequilibrium instabilities also in other sys- 
tems. For instance, a similar behavior has been observed 
in numerical studies of gauge field evolution for plasma 
instabilities in the context of relativistic heavy-ion colli- 
sions [2], where only the v dependence was studied. If we 
also restrict to the case of vanishing transverse momen- 
tum, the condition ( 23 1 yields the onset of exponential 
growth at times 



T Start 
To 



2/3 



for 



PT^O. (26) 



For these modes the subsequent growth continues as long 
as the linearized description remains valid. However, by 
investigating the entire momentum space one finds that 
there are modes which behave quite differently. Modes 
with vanishing longitudinal momentum exhibit a reso- 
nant amplification only until the time 



^TEnd 



To 



2/3 



for 



= 0. 



(27) 



corresponding conformal times by (10 1. This means that 



After this time these modes freeze out and exponential 
growth stops, while they start showing stable oscillatory 
behavior. 

The generic situation, where neither longitudinal nor 
transverse momentum vanish, is discussed in detail in 
appendix C and we refer here only to some characteris- 
tic results. By searching for real positive solutions for 9 
of the resonance criterion ( 23 ) as relevant for the initial 
value problem, we find that these solutions exist for all 
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FIG. 3. (color online) Time evolution of the equal-time trans- 
verse fluctuations F±{t, t,pt, v) from the linearized evolution 
equations for pt = and different rapidity wave numbers v. 
The modes are averaged over one period of oscillation of the 
macroscopic field. The time axis is scaled as (t/to)"^^^ . The 
vertical dashed grey li nes represent the set-in time of the in- 
stability according to (261. The full grey line corresponds to 



FIG. 4. (color online) Time evolution of the transverse fluc- 
tuations F±{t,t,pt,i^) from the linearized evolution equa- 
tions for 1/ — and different transverse momenta pr- The 
modes are averaged over one period of oscillation of the 
macroscopic field. The time axis is scaled as (t/tq)'^^^. The 
vertical grey dashed lines correspond to the estimates for the 
freeze-out times of the instability according to ( 27 1 . The full 



the maximum growth rate 3ffoT"o7o. 



grey line corresponds to the maximum growth rate 3(Toro7o. 



modes satisfying the time-independent condition 

1 



< 



54 



(28) 



Accordingly, all transverse modes {F±) satisfying the 
condition ( 28 ) experience exponential amplification for 



a certain period of time. The set-in and freeze-out times 
for these modes correspond to times where the relation 
( 23 1 is taken as an equality. Calling ^ = (54 Pt'^'^) / {'^o'''o) 
these times are given by 



1 + 2 sin I - arctan 



''"start 


^2/3 








QPt 


/ T'End^ 


^2/3 




V To y 




QPt 




1 + 2 cos I - arctan 



where f < 1 and, of course, rgtart < '''End for unstable 
modes according to (28). 



The phenomena of delayed set-in and freeze-out can be 
observed in our numerical studies of the linearized evolu- 
tion equations. In Figs. [3] and |4] we present the numerical 
solution of the linearized evolution equations for trans- 
verse modes F±{t,t,pt, v) with different transverse and 
longitudinal momenta. From Fig. [4] one observes that 
modes which are dominated by their transverse momen- 
tum exhibit an amplification at early times. This ampli- 
fication stops when they shift out of the resonance band. 
We find that this behavior is indeed very well described 
by our analytic estimate (27). In contrast, modes which 



are dominated by their rapidity wave number v, as shown 
in Fig. |3j exhibit an amplification at later times, when 



they have shifted inside the resonance band. Here the 
set-in time of the instability is well described by ( 26 ) . In 
both cases we find that the growth is exponential in con- 
formal time 9 oc t'^^'^ but with a time-dependent growth- 
rate 7(p(6')) < 7o. 

For instance, if we consider modes with small trans- 
verse momentum, such as pt — 0.2i7o and = as de- 
picted in Fig. |4j one observes that the growth rate is small 
initially, then closely approaches the maximum value of 
3o'o7'o7o and subsequently decreases until freeze-out oc- 
curs. The maximum growth ratc[^ 3fTo''o7o is realized 
only when j?{9) — pg, whe re jP p is the maximally am 
plified momentum given in (1251). Hence \i p^{9q) < 
initially, the modes experience a smaller growth rate 
at early times. For the considered mode p'^{9) ^ 9p\ 
and, therefore, it increases with (conformal) time. Hence 
l{p{^)) increases until it reaches the maximum rate when 
p{9) = po and subsequently decreases until the mode 
shifts out of the resonance band and the freeze-out oc- 
curs. This can be thought of as traversing the growth- 
rate 7(p) - as shown graphically in Fig. [12] in appendix 
B - from smaller to higher momenta (left to right in 
Fig. 12) in time. In contrast, modes with = and 
non- vanishing i/, as shown in Fig.[3j are characterized by 
a p{9) which decreases with time. Therefore, one finds 
that the largest growth rate is realized shortly after set- 
in of the instability and decreases towards later times. 



pI 



^ The maximum growth rate for equal-time correlation functions 
F{d, 9, Pt, v) is given by 270. In terms of the proper-time variable 
T this translates into 3a"of"070- 
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This corresponds to traversing the growth-rate 7(p) in 
the opposite direction, i.e. from higher to smaller mo- 
menta (right to left in Fig. 12 1 in time. 



rapidity 



III. 



NONLINEAR AMPLIFICATION OF 
INSTABILITIES 



A. Nonlinear corrections and power counting 

How to take into account nonlinear corrections in 
scalar quantum field theories is well known and can be 
based efficiently on the two-particle irreducible (2PI) ef- 
fective action , which circumvents problems of secular 
time evolutions encountered in non-resummed (IPI) ap- 
proximation schemes [5]. For the iV-component scalar 
quantum field theory a non-perturbative description can 
be based on the 1/A^-expansion of the 2PI effective ac- 
tion to next-to-leading order (NLO), which is explained 
in Refs. [T^]. The evolution equations have been pre- 
sented for general metric in Refs. [13l[14]. Gauge theo- 
ries with expansion have also been considered in Ref. [15] 
based on a classification of the 2PI effective action in the 
number of loops. 

Here we discuss those aspects of the time evolution, 
which are relevant for an understanding of the nonlin- 
ear amplification of instabilities. Nonlinear corrections 
will lead to strongly enhanced 'secondary' growth rates, 
which are multiples of the initial 'primary' growth rates 
observed from the linear regime in Sec. |Tl] Remarkably, 
this turns out to be very similar to the non-expanding 
case, which has been discussed extensively in the past 
for scalars O [9] as well as pure gauge theories [3| . 

The linear regime discussed above is described in terms 
of the macroscopic or background field (/), defined in Q, 
and the fluctuation or statistical two-point function F 
as given by Q. Together, the corresponding evolution 
equations ^ and ( 19 1 form a closed set of equations for 
the linear regime. An important new ingredient going 
beyond the linear regime will be the appearance of the 
time-dependent spectral function as an additional lin- 
early independent correlation function in the evolution 
equations. In general, for the considered scalar field the- 
ory there are two linearly independent two-point func- 
tions, which may be associated to the anti-commutator 
of two fields (F) and the commutator expectation value: 



Pabir, t',xt - x'rp,r]- r]') 

= i( 4>a{T,XT,v),M'r'i^'T,'n') 



(30) 



PabiT,T',PT,v)\T=T' 0, 

dTPah{T,T' ,pT,iy)\r=T' — 

T 

drdr'PabiT,T',PT,v)\r=T' ==0. 



(31) 



The r-dependence of the commutator between the field 
and its conjugate momentum enters via the metric tensor, 
y/—g = T [131 m] . Since these relations are valid at all 
times, they also fix the initial conditions for the evolution 
of pab{T-,T' ,pt,v). Again, using 0{N) symmetry we can 
write pab = diag(/9|[,/9_L, ...,pj_). 

Going beyond the linear regime using, e.g., the 2PI 
l/A^-expansion to NLO the field 0, the fiuctuations i^j^ || 
and the spectral functions p_Lj| form a closed set of cou- 
pled evolution equations. The linear equations ( 19 1 are 



generalized to their non-linear form, which reads for the 
longitudinal fluctuations 



a^ + i9.+p|+^ + Mf + A^^(r) 



- f dT"T" I]f[(T, t",pt, v) F|| (r", r',pT, v) 

J To 

+ r dT'W''j:f{T,T'',PT,iy)pil{T'',T',pT,iy). (32) 
J Tn 



Here the effective mass term = M^{F±^\\) and the 
non-zero spectral and statistical parts of the self-energy 
Sy'^ = Ey'^(pj^ II , _Fl II , 0) make the evolutions non- 
linear]^ The explicit linear T"-term in the integrand 
stems from the determinant of the metric tensor for the 
comoving coordinates. The spectral functions obey a 
similar equation with the characteristic 'memory inte- 
grals' over time, which reads for the longitudinal compo- 
nents 



1 



d^ + -dr+pi^ 



M|f 



2^ 

2N 



P\\{t, t',pt,v) 



dr T I]|[(t,t ,pT,y)p\\{T ,T,pT,iy). (33) 



The equivalent equations for the transverse components 
F± and p± can also be o btained from the corresponding 
linearized equations ( 19 ) by replacing with an effec- 
tive mass term -M^(^j.,|t) and taking into account a non- 



zero right hand side. The latter is of the same form as in 
( 32 1 and ( 33 ) with all longitudinal components replaced 



by transverse ones. At NLO in the 2PI 1/iV expansion 
the remaining equation for the field </> can be written in 



Here p denotes the spectral function determined by the 
commutator [.,.]. Therefore, it encodes the equal-time 
field commutation relations, which read in Fourier-space 
with respect to the transverse spatial coordinates and 



* The spectral part, E^, can be related to the imaginary part and 
the statistical part, S^, to the real part of the self-energy for the 
considered theory [9]. 
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FIG. 5. (color online) Vertices in the presence of a macro- 
scopic field. Dashed (red) lines denote 'longitudinal' field 
components and solid (black) lines are associated to either 
longitudinal or 'transverse' components. The dot indicates a 
non-zero field expectation value. 



the form 



dr'r' I]||(t, t',pt = 0,u^ 0)U=o 0(t') , (34) 



where the spectral part of the self-energy is evaluated for 



zero field, i.e. (p_|__|| , F|_ y 



0) [Hlin]- Of course, at 



even higher order in 1/N also additional terms that are 
non- linear in the field (j) appear |12) . 

We will start by classifying the nonlinear corrections 
to the evolution equations entering via the effective mass 
or self-energy terms || , S'J^ y and S;^ || . In order to 
write down self-energies, it is important to note that in 
addition to the four-vertex proportional to A there is an 
effective three-vertex due to the presence of the macro- 
scopic field. This is visualized in the left panel of Fig. [5] 
where dashed (red) lines denote 'longitudinal' (a = 1) 
field components and solid (black) lines are associated to 
either longitudinal or 'transverse' (a = 2, . . . , TV) compo- 
nents. A dot indicates a non-zero field expectation value 
(^(t) as in 

A general non-linear contribution to the evolution 
equations will contain powers of A, of the 'propagators' 
Fl II and pj.,||, and of the field (p. Initially, the para- 
metric dependence of the field is ~ 1/A. Therefore, 
at not too late times a classification based on a small A 
has to take this into account. Most importantly, taking 
-FljI into account for the power counting is crucial, since 
F grows exponentially in time as a consequence of an in- 
stability. In contrast, the 'weight' of the spectral function 
Pj^ II remains parametrically of order one at all times as is 
encoded in the equal-time commutation relations (31 ) 



It is also important to note the fact that transverse fluc- 
tuations {Fj_) exhibit the dominant growth in the lin- 
ear regime. Consequently, contributions containing more 
transverse propagators {F±) can become important ear- 
lier than those diagrams containing longitudinal propaga- 
tors (-Fy) instead. For instance, an expression containing 
powers A"F™0^' with integers n, m and I may be ex- 
pected at not too late times to give sizable corrections to 
the linearized evolution equations once 
for typical momenta. Here n yields the suppression factor 
from the coupling constant, whereas m introduces the en- 
hancement due to large fluctuations for typical momenta 



and / due to a large macroscopic field. The power count- 
ing can become more involved as time proceeds, and it 
is remarkable that one can indeed identify a sequence 
of characteristic time scales with corresponding growth 
rates. 



B. Characteristic time scales and growth rates 

To start with a simple example, we consider first one- 
loop 'tadpole' corrections, which are obtained by 'closing' 
two longitudinal or two transverse 'legs' of the four- vertex 
in Fig. [5] on the right. We emphasize already here that 
there are other corrections which will be of relevance be- 
fore tadpoles come into play for not too large N. How- 
ever, it will be convenient to express time scales in terms 
of the characteristic time when tadpoles become relevant, 
since this turns out to coincide with the time when an 
infinite scries of corrections become sizable and all expo- 
nential growth of fluctuations stops. Tadpole corrections 
are mass-like and their contribution reads 

M2(r) = m2 + T||(r) + (Ar + l)T_L(T), 

M^{t) = m^ + 3T||(t) + {N - 1)Tj_{t) (35) 

for the nonlinear evolution equation of the transverse 
(Fl) and the longitudinal fluctuations (Fy), respectively. 
The one-loop tadpole integrals read 



T±,||(r) = 



2^ 



cPpTdv 
(27r)3 



F±niT,T,pT,iy) , (36) 



where A denotes some suitable regularization that enters 
the renormalization procedure, which has not to be spec- 
ified for the current purpose. While initially the tadpole 



contributions ( 36 ) are suppressed by the coupling con- 



stant A, they become sizable at later times once the fluc- 
tuations grow parametrically to Fj^ || (t, T,p7-, i/) ^ 1/A 
for typical momenta pt and v. More precisely, for the 
massless case (m^ — 0) the size of a tadpole contribution 
in ( 35 ) should be compared to that of the macroscopic 



field-squared term A<^ (t)/67V in (34), which on aver 



age is (7^/2 (t/tq)"^^'^ at early times. Wc will denote 
the time when both become of the same order of magni- 
tude as Tnonpcrt- Using that the tadpole will quickly be 
dominated by unstable modes entering the integral with 
characteristic primary growth-rate 7(^7, 1^,0) < 70, such 
that 



F{e,e',pT,v) 



F{eo,eo,PT,j^) 

one obtains the estimate 
1 



< exp [70 (9 + 9'- 29o) 



"^nonpcrt ^ '^0 



1 



-In 



3(7oro7o V2(A^ + 1)T±(to) 



(37) 



3/2 



(38) 



Here we used the important fact that transverse fluc- 
tuations exhibit the dominant growth in the linearized 
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FIG. 6. (color online) One-loop contribution to the longi- 
tudinal component of the self-energy. Transverse two-point 
functions are denoted by dotted (blue) lines. 



evolution equations. It is noteworthy that Tnonpcrt is 
rather sensitive to the inverse of the primary growth-rate, 
whereas the coupling constant and the size of the initial 
fluctuations only enter logarithmically through rL(To) 
according to (36). In the weak coupling limit (38) re- 
duces to 



' nonpcrt 



> 



To 



-ln( - 
3(Toro7o V A 



3/2 



(39) 



at leading logarithmic accuracy. 

We strongly emphasize that when fluctuations have 
grown parametrically to 0{\/\) at Tnonpcrt, it is not only 
a one-loop tadpole that becomes of order one. Since there 
is an inflnite series of contributions including arbitrarily 
high loop-orders that become sizable, a non-perturbative 
approach such as the 2PI 1/7V expansion or classical- 
statistical simulation methods have to be applied. Before 
addressing the non-perturbative regime in Sec. |IV| below, 
where no longer an exponential growth will be observed, 
we first consider earlier times than Tnonpcrt to discuss the 
nonlinear amplification of instabilities. More precisely, 
we will focus here on a subset of diagrams of the 2PI 
1/N expansion to NLO [HI, which are relevant for times 

T ^ Tnonpcrt • 

In general, the smaller the above introduced param- 
eter (n — l)/ra for a specific self-energy correction, the 
earlier it may be expected to play a sizable role during 
the nonequilibrium time evolution. For the tadpoles this 
parameter is one. However, there is another one-loop self- 
energy correction for which this parameter is 1/2, which 
solely contributes to the longitudinal components of the 
self-energy, as indicated by the red dashed (amputated) 
legs. The respective one- loop contribution to the longi- 
tudinal component of the self-energy is displayed graph- 
ically in Fig. [6j Here we denote transverse two-point 
functions by dotted (blue) lines. It is important to note 
that in this power counting scheme all other diagrams are 
suppressed by at least a fractional power of the coupling 
constant A. Hence there exists a kinematic window where 
the only relevant self-energy correction originates from 
the diagram shown in Fig. [6j In the following, we will 
see that this leads already to a non-linear amplification of 
the primary instability for longitudinal fiuctuations. In 
particular, this correction has to be taken into account 
when considering further corrections that become rele- 
vant at later times. 

The spectral and statistical self-energies associated to 



Fig. |6] read 

p „ (one-loop) 4A(A^ - 1) „ 

2j||(t, T tPt^v) — — — a(T)a[T ) 



(27r)3 



6iV 

P_l(t, t",pt -qr.v- i^q)Fj_{T, t", qt, Vq) 



(one-loop) 2A(iV- 1) 



dPqTdvq 



(27r)3 



F^{t, t",Pt -qr^V - Vq)F^{T, t" ,qT,Vq) 

(40) 



-^PA^{t, t",Pt -qT.V- l^q)p±{T, t" ,qT,Vq) 

These enter the memory integrals on the right-hand-side 
of the evolution-equations ( 32 ) , ( 33 ) and ( 34 ) . We note 
that the (pj^pj^)-term in the integrand for Ey is a genuine 
quantum correction, which would be absent in a classical- 
statistical description 9 . However, since F^F±^ ^ 
once nonlinear corrections become sizable, one can ne- 
glect the quantum part to very good accuracy. We will 
exploit this fact further in Sec. |IV| on classical-statistical 
simulations below. 

In order to make analytical progress, one can exploit 
the fact that the dominant contribution to the mem- 
ory integrals originates from late times when fiuctuations 
have become exponentially large [5] . Instead of consider- 
ing integrals from tq to t and t', respectively, we consider 
the integrals only over some suitable, small interval A. 
This will be sufficient to obtain characteristic time scales 
to leading logarithmic accuracy. It allows us to expand 
the integrand around the times of interest, where at lead- 
ing order one finds 



Fir, t",pt,v) ~ F{t, t,pt,v) , 

I N ^ t-t" t- 
P(t,t ,Pt,v) ^ ^ 7. 



(41) 
(42) 



Here we used the equal-time commutation relations (31 ) 



to expand the spectral function. With these approxima- 
tions one can explicitly evaluate the right-hand-side of 



the evolution equation ( 32 ) as 
1 



RHS 



3N 



CpqxdVq 



xF||(t,t',p^,.) + A2 



6iV 



(27r)3 



F±{t, T, qT,l'q) 

(43) 



d^qxdvq 
(27r)3 



F±{t, t',Pt - qr.V ~ Vq)F^{T, t' .qr^Vq) ■ 



The first term in the above sum acts as a momentum- 
independent mass term similar to the tadpole term on the 



LHS of the evolution equation (32 ). For it to be relevant 



it requires Fi_ ^ 0{l/\) as discussed before and, hence, 
one expects a relevant contribution staring at Tnonpcrt- 
In contrast, the second term appearing in (43) acts as a 



source term for the evolution of the longitudinal fluctu- 
ations i^ij . The momentum dependence is given by the 
convolution of the transverse fluctuations F±^ with itself. 
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FIG. 7. (color online) One- and two-loop contributions to 
the transverse components of the self-energies. 



For instance, if F±(T,T,qT,Vq) is peaked around some 
momenta po and ±i^0: the source term has its dominant 
contributions around px = {0, 2po} and i' = {0, ±2i'o}- 

In the presence of this strong source term, one ex- 
pects F|| to follow the source ~ A-Fj. In this way the 
primary growth of the transverse fluctuations F± leads 
after some delay-time Tgource to a secondary stage of 
growth, where the longitudinal modes Fy exhibit ampli- 
fied growth-rates. In particular for the longitudinal zero- 
mode, i.e. F'|j(r, r',pT = 0, 1/ = 0), the secondary growth- 
rate can be twice as large as the primary growth-rate 
observed for F±. It is important to note here that the 
amplification occurs before tadpole corrections become 
important. This is due to the fact that the source term 
is of order XF]_ which requires F± ^ C(l/Ai/2) for typ- 
ical momenta to yield relevant contributions. With the 
same approximations entering the estimate for Tnonpcrt 
one can estimate the associated time-scale as 



^source ^ ^0 



3/2 



(44) 



Again, in the weak coupling li mit the dominant contri- 
bution arises from ln(A^^) and (44) reduces to 



(-''^1) ■''nonpcrt 
- 23/2 



(45) 



which is smaller than the time-scale on which screen- 
ing effects due to effective mass terms become relevant. 
Hence, there is a period of time when one expects growth 
of the longitudinal modes due to a nonlinear amplifica- 
tion of the primary instability. This happens for a bound 
momentum region which is entirely determined by the 
spectral shape of the primary instability. We empha- 
size that the larger couplings and/or initial fluctuations 
the earlier this non-linear amplification of the instability 
happens. 

We have seen that as a consequence of nonlinear am- 
plifications characteristic longitudinal fluctuations can be 
expected to become 0{XF^) around the time Tgourcc- 
This modifles the power counting for subsequent times, 
since in addition to parametrically large Fj_ also para- 
metrically large -Fy enter loop corrections. In particu- 
lar, this will soon after lead to sizable nonlinear contri- 
butions to the evolution equation for transverse fluctu- 
ations F±. The relevant diagrams contributing to the 
transverse components of the self-energies are displayed 
in Fig. [Tj Assuming F|| ^ ©(AFj) for typical momenta. 
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FIG. 8. (color online) Time evolution of transverse fluctua- 
tions with pt = 0, 4, 8, 12, 16, 20 pmin and 1^ = averaged over 
one period of oscillation of the macroscopic field. 



both depicted diagrams are ©(A^Fj). Thus one ex- 
pects a relevant contribution as soon as Fj_ ^ l/A^/^. 
The momentum-dependence of the diagrammatic con- 
tribution leads to an extension of the amplified region 
to higher momenta, where again multiples of the pri- 
mary growth-rate appear. It is important to realize that 
this amplification repeats itself, i.e. the newly amplified 
modes together with the primarily amplified ones act as 
a source for other modes. In this way the instability 
propagates to higher and higher momentum modes for 
both longitudinal and transverse fiuctuations. This can 
be nicely observed from our numerical simulations, which 
are presented in the next section. 



IV. CLASSICAL-STATISTICAL LATTICE 
SIMULATIONS 

A. Comparison with analytic results 

As the system evolves in time, occupation numbers 
may grow until the fluctuations Fj^ || become as large as 
0{1/X) around the time T„onpcrt given by (38). Hence, 
there are sizable contributions to the dynamics originat- 
ing from all loop orders and a non-perturbative descrip- 
tion is necessary. The growth saturates at this time and 
a comparatively slow, quasistationary evolution sets in. 
This has been discussed in great detail for non-expanding 
systems ^ . In that context it has also been shown that 
classical-statistical approximations can reproduce very 
well the results from the non-perturbative 1/N expan- 
sion to NLO of the quantum 2P1 effective action at not 
too late time, i.e. before the approach to quantum ther- 
mal equilibrium sets in [16H18) . 

In this section we present the results obtained from 
classical-statistical lattice simulations for the longitudi- 
nally expanding case discussed above. These are de- 
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FIG. 9. (color online) Time evolution of longitudinal fluctu- 
ations with pt ~ 4, 8, 12, 16, 20 Pmin and u — averaged over 
one period of oscillation of the macroscopic field. 



FIG. 10. (color online) Time evolution of transverse fluc- 
tuations with Pt = 4 pmin and = 0, 10, 20, 30, 40, 50 Umin 
(top to bottom) averaged over one period of oscillation of the 
macroscopic fleld. 



scribed by numerically solving the classical field equa- 
tions of motion, obtained from the stationarity of the 
action ([T]), and Monte Carlo sampling of initial condi- 
tions such that the initial classical averages agree with the 
quantum initial correlation functions. The parameters of 
these simulations are chosen as = 0, a^TQ = 5, N = A 
and A = lO^"*. We discretize the evolution equation for 
the inhomogeneous classical field on a three-dimensional 
lattice in transverse coordinates and rapidity with grid 
size Nf X iV^. The lattice spacing is chosen such that 
the lattice ultraviolet cut-offs ~ l/ot and ~ l/(Ta^) are 
above all physical scales. Of course, for the expand- 
ing system this condition is time-dependent through the 
time-dependence of the lattice cut-off due to red-shift as 
well as the time-dependence of physical scales due to dilu- 
tion. Numerical simulations at late times are, therefore, 
computationally hard to perform and we will focus the 
discussion on the physics at sufficiently early times. If 
not stated otherwise we will use Nt = N,^ = 128 with the 
discretization aoat = 0.5 and = 0.1. For convenience 
we express the transverse and longitudinal momenta in 
terms of the lattice momenta Pmin = 2y/2/{Ntat) and 
i^min = 2.0/(iV^a^), respectively. 

The time-evolution of transverse and longitudinal fiuc- 
tuations for different momenta is presented in Figs.|8]and 
[9] for different transverse momenta px and rapidity wave 
number = 0. Fig. [lO] shows the transverse fluctuations 
for modes with different rapidity wave numbers v and 
transverse momentum = 4 Pmin for which the leading 
primary instability occurs. From Fig. [8] one observes that 
the primary instability occurs for a bound momentum 
region. The phenomenon of linear freeze-out is clearly 
visible. The modes with px — S Pmin and 12 Pmin, for in- 
stance, exhibit exponential growth at early times but de- 
couple from the instability shortly. The functional form 
of the primary instability is well described by an expo- 
nential in (t/tq)'^/^ with maximum growth-rate 3cro''o7o- 



At later times, when the primary instability has been 
operative long enough to produce large transverse fluctu- 
ations, secondary instabilities set in for the longitudinal 
modes, as can be observed from Fig. [9j These secondary 
instabilities exhibit growth rates up to SfJoToTo a-s dis- 
cussed in Sec. |III[ The onset of these secondary insta- 
bilities is limited to a momentum region pt ^ 8 Pmin 
as they originate from the two-loop diagram shown in 
Fig. [7| This first non-linear amplification of the primary 
instability happens at times t/tq ~ 15. At later times 
t/tq w 20 the system exhibits collective amplification of 
the primary instability and fiuctuations begin to grow in 
a wide momentum range. Also modes that previously ex- 
hibited freeze-out show a second period of growth while 
modes with smaller primary growth-rates exhibit a sig- 
nificant speed-up. One can also observe from Figs. [8] and 
|9]how this sets off an avalanche of instabilities propagat- 
ing to higher momenta. At this point of the evolution, 
modes with high rapidity wave number v set in earlier 
as suggested by the linearized evolution equations. From 
Fig. [To] one observes that once amplification by nonlinear 
corrections sets in, modes in a large momentum region 
show significant growth within a very short time. This 
happens because high longitudinal momenta no longer 
have to wait to become unstable, but are instead sub- 
ject to the nonlinear amplification. This is indeed very 
similar to what is observed in numerical simulations of 
pure gauge theory with 'Glasma' initial conditions f2]. 
Ultimately the growth of the instability saturates when 
the fluctuations become of the order of the inverse cou- 
pling 0{1/X). The rapid dynamics of instabilities is then 
followed by a regime of comparatively slow and smooth 
evolution, where one expects transport approaches to be 
applicable. 
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B. Energy transfer to particles and EOS 



^particles''^ 



Finally one is interested in the applicability of hydro- 
dynamics to provide a simple description of the system 
under consideration. A necessary condition for this to 
be the case is the establishment of an equation of state 
relating pressure and energy density. To investigate this 
in our numerical simulations we consider the different 
components of the stress energy tensor 

T^^(x) = WM^)) id''^a{x))-crCMx)) . (46) 

The quantities of interest here are the energy density e 
as well as the longitudinal and transverse pressures, Pl 
and Pt, defined byj^ 

e{T) ^ ^ J d'^XTdr]T''^{T,XT,v) , 
Pt{t) ^ y J ^ ^Tdi-j , 

PL{r) = yj d^XTd7jT'^'^{T,XT,v) ■ (47) 

To further investigate the energy transfer from the 
macroscopic field to the fluctuations, we define the en- 
ergy content of the macroscopic field by replacing the 
fluctuating field in (46) by its expectation value, i.e. 



T^-ix) = (d'^M^)) id^M^)) - g^'^mi^) ■ (48) 

The energy contained in the fluctuations is then simply 
the difference between the two, i.e. 



T'"' (t) 

^ particles V-^^ 



-T^''{x)-Ti;%x) 



(49) 



In the upper part of Fig. [TTJwe show the energy fractions 
contained in the macroscopic field and the produced par- 
ticles as a function of time. One observes that the en- 
ergy contained in the produced particles becomes signifi- 
cant around t/tq ~ 20, which coincides with the onset of 
nonlinear amplification in a broad momentum range dis- 
cussed above. While the large population of soft particles 
is created at early times through the primary instability 
as sugg ested by Figs. [8] and ^ the energy transfer to 
the particles is to be understood as a consequences of 
hard particles being affected by secondary instabilities 
through collective amplification because of the enhanced 
phase-space and larger mode energy for hard particles |^ 
The ratios of longitudinal and transverse pressure over 
energy density, i.e. the dynamical equation of state, are 
presented in the lower panel of Fig. [TT] as a function of 



Since we assumed a boost-invariant system the local restframe 
is always given by uf^{x) = (cosh{?;), 0, 0, sinh(7y)) Therefore, one 
obtains energy-density and pressure immediately in terms of the 
components of the stress-energy tensor in comoving coordinates. 
In a pertrubative treatment one has T^^^{X) = 

J -j^^^p'^p'^ F{X,p), where F(X,p) is the Wigner transform of 

the statistical two point correlation function [9]. 




10 20 30 40 50 60 70 80 90 100 

X/Xn 




10 20 30 40 50 60 70 80 90 100 

X/Xn 



FIG. 11. (color online) (top) Fraction of total energy con- 
tained in the macroscopic field and in the produced particles. 
Significant energy transfer to the produced particles occurs 
when high momentum modes carrying more energy are af- 
fected by secondary instabilities due to nonlinear amplifica- 
tion in a broad momentum range, (bottom) Dynamical equa- 
tion of state for transverse and longitudinal pressure. 



proper-time. Initially this observable is dominated by 
the energy contained in the macroscopic field. Therefore 
one observes large oscillations around the average value of 
1/3. At later times when energy is carried predominantly 
by the produced particles the oscillations are damped and 
one observes a smooth behavior. In our simulations we 
find that the system exhibits a remaining anisotropy as 
longitudinal and transverse pressure do not equilibrate to 
1/3. It is important to realize here that isotropy can only 
persist when the evolution of the system is dominated by 
particle interactions in contrast to the free-streaming be- 
havior. The latter drives the system away from isotropy 
due to the red-shift of longitudinal modes. These quan- 
titative aspects depend, of course, on the specific model 
and can be very different in the relevant case of QCD. 
At late times one should also employ the quantum evolu- 
tion based on the 2PI 1/A^ expansion to NLO instead of 
classical-statistical simulations, since deviations can be 
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expected for a dilute system. Unfortunately, there is no 
corresponding nonperturbative expansion so far that is 
manageable in QCD and this topic is beyond the scope 
of the present work. 



V. CONCLUSION AND OUTLOOK 

We have shown for the example of A^-component scalar 
field theories that nonequilibrium instabilities can lead 
to important nonlinear phenomena due to quantum cor- 
rections for longitudinally expanding systems. At early 
times fluctuations are amplified due to primary insta- 
bilities, which are described by the linearized evolution 
equations. When these become large enough to induce 
nonlinear corrections, the instability propagates towards 
higher momentum modes and secondary growth-rates 
much higher than the primary ones can be observed. This 
way significant energy densities carried by the produced 
particles can be achieved, which has important conse- 
quences for the question of isotropization and the estab- 
lishment of an equation of state. 

The presence of secondary instabilities is rather insen- 
sitive to the detailed mechanism for primary growth of 
fiuctuations. Qualitatively, very similar phenomena can 
be expected for nonabelian gauge theories as confirmed 
already for the non-expanding case [3l 110) . This has in 
principle been observed also in classical-statistical sim- 
ulations of the Glasma [H [S] and corresponding 2PI ef- 
fective action techniques |15| may provide a systematic 
framework to study these effects. Even for the physics of 
primary instabilities with anisotropic Bjorken expansion 
many qualitative aspects, such as the delayed set-in of 
high v modes observed in classical-statistical lattice sim- 
ulations in SU(2) Yang-Mills theory [2", can also be found 
in the scalar theory and explained analytically within our 
model. 

The relevant time-scale for nonlinear effects is governed 
by the growth rates of the primary instability and to a 
weaker extent by the logarithmic dependence on the cou- 
pling constant and the spectrum of initial fluctuations. 
The evolution of the Glasma in the context of relativistic 
heavy ion-collisions may be strongly modified once real- 
istic initial fluctuations are taken into account [TH] ■ Here 
one expects the instability, which affects initially only 
small longitudinal momentum modes, to extend quickly 
to higher momenta by the nonlinear amplification mech- 
anism if their nonzero initial amplitude is properly taken 
into account. 

In this work, we have not discussed the possible 
approach to a turbulent scaling regime once nonlinear 
corrections become of order one. It has been argued in 
the context of QCD that also a Bose condensate may 
form later on |2D]- For the non-expanding case it has 
recently been shown for scalar theories that indeed a 
Bose condensate occurs as a consequence of an inverse 
particle cascade with a universal power-law spectrum 
pn [55] . This particle transport towards low momenta is 



part of a dual cascade, in which energy is also transfered 
by weak wave turbulence towards higher momenta. 
Similar phenomena have also been studied for cold 
atoms [51] [53] . Classical-statistical simulations including 
the relevant physics of longitudinal expansion [51 [S] 
should be the appropriate theoretical approach to clarify 
these questions in the context of heavy-ion collisions at 
sufficiently high energies. 

This work is supported in part by the BMBF grant 
06DA9018 and by EMMI. 



APPENDIX A: CONFORMAL TIME IN 
ANISOTROPICALLY EXPANDING GEOMETRY 

Here we introduce the conformal time on the level of 
the classical action. This has the advantage that the scale 
factor appearing in the 2PI evolution equations (32 1-(34) 



can be absorbed into a redefinition of the time-variable 
and a rescaling of field variables and propagators. In con- 
trast to the main text we do not make the conformal time 
variable dimensionless. The classical action in comoving 
coordinates is given by ([T]) 



d'xT drj dT air) -.g^'^(x) (d^ipa) (d^ifa) 



1 



X 



(50) 



with the metric g^'^{x) — diag(l, — 1, — 1, — r^) and the 
scale factor a(r) which in the case of a one-dimensional 
Bjorken expansion is given by the determinant of the 
metric a(T) = —g^u{x) = t. In isotropically expanding 
backgrounds it is a frequently used concept to introduce 
conformal time variables to solve the dynamics in a quasi- 
static framework. Here we employ an approach that is 
very similar in spirit. By introducing a new time variable 
and rescaling the field variables according to 



dT 



,1/3 



{T)de 



,1/3 



(51) 



where 9 will be referred to as conformal time, the action 
can be rewritten as 



d'xT dfl d0 { ^r^i^) (d^^^a) (d^ifia) 



(52) 



where the determinant of the metric drops out 
and first order derivatives are absorbed into a re- 
deflnition of the effective mass through integration 
by parts. Here the metric is given by g^^ = 
diag(l, -a-2/3(-r)^ _q-2/3^^)^ _a4/3^^))_ j.^^ ^^^.^^^ 

appearing in ( |52[ ) depends explicitly on conformal time 
and is given by 
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where primes denote derivatives with respect to confor- 
mal time. The last two terms originate from interchang- 
ing time-derivatives with the scale factor 



n2 l,nl~J la'W- 



rfr a(T) (driPa) = I de[tf^- 



,(54) 



and integrating first derivatives by parts according to 

'a"{0) (a'{ey 



{6) V a{e) 



The scale factor and its derivatives with respect to con- 
formal time are explicitly given by 




1/2 



-1/2 



(56) 



The major advantage of the formulation in terms of con- 
formal variables is the fact that the metric determinant 
drops out of the action. Therefore the evolution equa- 
tions derived from the action ( 52 ) do not contain first 



order time derivatives. However the mass term m?{0) 
depends explicitly on time and dominates on large time- 
scales, indicating the freeze-out, unless m? = 0. In con- 
trast for the massless case one finds that Tn^{6) decreases 
as a power of on large time-scales. Hence the confor- 
mal time formulation is particularly well suited for this 
situation, which is relevant also for gauge theories. 



APPENDIX B: SOLUTION OF THE LAME 
EQUATION 

Here we give some details for the solution of the 
Lame equation as it appears in parametric resonance in 
Minkowski space-time needed for the main text. The 
discussion is presented for the massless case = and 
closely follows in part Ref. [Mj . In this situation the lin- 
earized evolution equation for the fluctuations for trans- 
verse modes is given by 



di 



p 



F{e,e',p)^o. (57) 



By introducing the mode functions fp{9) according to 

Fi0,9',p) = I [U{e)f;{0') + f;ie)u{9')] (58) 



the partial differential equation ( |57[ ) can be transformed 
into a set of two independent ordinary differential equa- 
tions for evolution in and 9' respectively. By use of the 
identities (23 



cn'^{9,a) = 1 - STi'^{9,a) , 



sn2(0,a) 



a sn"^ {9 + iK' {a), a) ' 



(59) 
(60) 



A 


ei 


62 


63 


52 


93 


1 


1/2 





-1/2 


1 






TABLE 1. Parameters of Weierstrass Elliptic functions for 
a = 1/2 



where K'{a) = K{l—a) and K{a) is the complete elliptic 
integral of the first kind, ([57]) can be represented in terms 
of Weierstrass functions p0) by virtue of 



p(^ ; 32,53) 63 + 



ei - 63 



(61) 



sn2(yei - 636', a) 

The freedom of choice of the prefactor of p(z ;52,53) 
allows us to write the evolution equation as 

[d^e+P^~MO + ^K'{a); 92, 93)] fp{9) = , (62) 

where the parameters for a = 1/2 are summarized in 
Tab. HI By expressing the momentum implicitly in 
terms of the Weierstrass function p(z) with 



p{z ■,92,93) = -P^ 



(63) 



( 62 1 reduces to the well known Lame equation in the 
Weierstrass form. The independent solutions Up{9) 
and Up{—9) are given in terms of Weierstrass functions 
C{S;92,93) and a{9;92,g3) and read 24, 26J 

rr (n. ^^ec(.) '^{0 + ^ + ^K'(a)) a{^K'{a)) 

' " a{9 + iK'{a)) a{z + iK'{a)) ' ^ ' 

To further discuss the properties of these solutions for 
different momenta p one needs to investigate the map- 



ping (63). The Weierstrass elliptic function p{z) takes 
real values on the sides of the fundamental rectangle. 
Furthermore as > we can restrict to negative values 
of p{z). The solution to (63) can be parametrized as 

z^iP forp2>l/2 (65) 

z = iK' {1/2) + /3 for < p2 < 1/2 (66) 

where < /? < ^(1/2). The growth-rates can be ob- 
tained from a Floquet analysis, i.e. by investigating 



Up{9 + 2K{l/2)) = e2^(i/2>^(P) Up{9) . 



(67) 



From (64) and using properties of the Weierstrass func- 
tions |25| one finds 



F{p) = 



1 



K{l/2) 



{zaK{l/2))-K{l/2)az)) . (68) 



If this has a non-vanishing real part then the solution 
exhibits exponential growth and belongs to the resonance 
band. This is the case for all modes with p'^ < 1/2. If 
in contrast pP' > 1/2, the Floquet index F{p) is purely 
imaginary. Then the mode exhibits stable oscillations 
and therefore belongs to the stable band. In particular, 
the stable modes do not exhibit growth and are thus not 
relevant relevant for our discussion of nonlinear effects of 
the evolution. We therefore concentrate on the properties 
of the resonance band. 
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B.l Resonance Band 

To further study the resonance band it is useful to 
invert the defining equation 



z{p) = iK' 



- m ■ 



(69) 



By use of the identities (59)-(61) and performing the 



steps (59)-(63) 'backwards' we find 
, 1 



P 



(70) 



for resonant modes satisfying the inequality ( 66 ) . To 



obtain a more intuitive expression for the growth-rates 
the Floquet indices F (p) can be expressed as 



Z{u; a) 



4 {2K{a) ) 



2K{a).o ' 



(71) 



where ^j{u), j — 1, 2, 3, 4 are the Jacobi theta func- 
tions. By use of the relations between Weierstrass zeta 
and Jacobi theta functions 



C(^) 

and using 



K(a) 2K(a)M ( 

^1 \2K(a) 



(72) 



al ( , 7T K'{a)\ 



it is straightforward to obtain the final result 



F{p) 



2K{l/2] 



- ZiPip): 1/2) 



(73) 



(74) 



The Jacobi zeta function Z(/3; 1/2) is purely real and 
therefore yields the growth-rate 

7(p) = Z(/3(p); 1/2), /3(p) = cn-i(yV). (75) 

The growth rate as a function of momentum is shown 
in the left panel of Fig. 12 The imaginary part in (74) 



is precisely the oscillation frequency of the macroscopic 
field. It is useful to apply approximations to the growth 
rate in ( 75 1 in order to obtain the explicit momentum 



dependence. To simplify the expression ( 75 1 we expand 
the inverse Jacobi cosine according to 



cn-i(z; a) = cos~^ (z) + 0{a) 



(76) 



In the next step we expand the Jacobi zeta function in a 
g-series 



Z(/3, 1/2) 



277 



if (1/2; 



/ 2nnl3 
\2K{l/2) 



(77) 
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FIG. 12. (color online) Properties of the solution of the 
Lame equation in Minkowski space-time, (left) Comparison 
of growth rates of fluctuations fro m th e exact solution of 1 68 1 
(black) with the approximation of (781 (red dashed), (right) 
Oscillation frequency obtained from the exact solution ( |68[ ) of 
the Lame equation (black). Also shown is the free field limit 
(red dashed) and the oscillation frequency of the macroscopic 
field (blue dashed). 



where q = e~'^ . By keeping the leading term in the above 
q-series and using trigonometric identities we obtain the 
final result 



lip) 



Aire 



K(l/2) 



V2p2 (1 _ 2p2) . 



(78) 



This is compared to the exact growth-rate in Fig. [12] 
While the overall behavior is described rather well, one 
observes that there are some deviations for the maximally 
amplified modes. For the latter one finds in the above 
approximation 



7o 



27r e-'^ 

mm 



for 



Po 



(79) 



16 



APPENDIX C: QUASI-STATIC 
APPROXIMATION 

In the quasi-static approximation discussed in 
Sec. |IIC| we assume that the time dependence of the 
momenta is slow on the time scale of one period of oscil- 
lation of the macroscopic field. The explicit time depen- 
dence of the momenta p{6) in the evolution equation for 
the fluctuations 

[d^g+f{9) + cn^{0~eo; i/2)]F^{e,e',PT,i^)^o 

(80) 

can therefore be neglected when inferring the momentum 
dependent growth-rate. In practice, this corresponds to 
replacing p p(jPt, 9) for all expressions presented in 
appendix B. We test this approximation below, where we 
present a comparison between numerical simulations of 



1 

0.8 
0.6 
0.4 - 
0.2 




( 80 1 and the quasi-static approximation. 



FIG. 13. (color online) Visualization of (84 1 for different 
values of pyi/^/((To (ctoto)^)- There is always one negative so- 
lution which is irrelevant for the initial value problem. If 
Pyi^^/(o-o(aoro)^) is smaller than the critical value of 1/54 
there are two positive solutions corresponding to the set-in 
and freeze-out times of the primary instabihty. 



C.l Set-in and freeze-out times in the quasi-static 
approximation 

By replacing p — p{pT,i',0) the condition for modes 
to be contained in the resonance band, becomes time 
dependent. For p^{e) = 2e/{3aoTo) p|/crg + 9^^/(4:9^) 
the condition then reads 



29 



Pt 



1 

3(cToro) a2 ' 49^-2 



(81) 



For the initial value problem we are only interested in 
the solutions in the forward light-cone, i.e. solutions for 
which 9 > 0. Before considering the generic case pr ^ 
and v we consider briefly the special cases where one 
of the two vanishes. For vanishing transverse momentum 
the condition is satisfied if 



(82) 



which suggests that exponential growth sets in with a 
delay. In contrast for vanishing ly growth is limited to 
the time when 



Sag 



(83) 



suggesting that exponential growth stops for later times. 
For the generic case (81 1 can be solved graphically. First 
we rescale 9' /9 ~ 2/(3ctoTo) p^/ctq yielding 



9 4 



1 



CTo(croTo) 



< 



1 



(84) 



This suggests that the existence of an unstable window 
only depends on the value of v'^p'^, whereas the pre- 
cise time also depends on p\. The existence of real 
solutions is visualized in Fig. |13[ where we show the 
LHS and the RHS of ^ for different values of i/^pf, . 



One observes that one solution is always real and neg- 
ative and therefore not relevant for the study of the 
initial value problem. The other solutions are real for 
v^Px < fo('''o'''o)^/54. In this case the solutions are given 

by 



^ Start 



9' 

'End 



_ 1 

~ 6 

4[i 



/3 = 



6 
54 



1 - e-"/Va - e 
+ l/a + a] , 

2/3 + (/? . 

4 2 

Pt^ 



1/3 



(85) 

(86) 
(87) 



The solution is real if and only if \a\ = 1. Indeed this 
is the case for resonant modes, i.e. modes with /3 < 1. 
We can exploit this fact to rewrite the solution for the 
starting and ending time as 



'start 



'End 



(1 -|-2sin(0/3- 7r/6)) 
(1 -|-2cos((?!)/3)) 



= 2 arctan 



1 



1 



(89) 
(90) 
(91) 



which is manifestly real. We note that (j) € (0, tt) for 
modes satisfying ( 84 ) and therefore 9'^^^^^^ < 9'^^^ as sug- 



gested by the naming. 



C.2 Tests of the approximation 

The quasi-static approximation yields growth rates 
which can be integrated in time numerically. In Fig. [14] 
this estimate is compared to the full time evolution of 
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FIG. 14. (color online) Time evolution of the transverse fluc- 
tuations Fx {t,t,pt,i^) averaged over one period of oscilla- 
tion of the macroscopic field. The black lines correspond to 
the numerical solution of the linearized evolution equations 
and red lines show the corresponding solutions in the quasi- 
static approximation where the growth rates have been in- 
tegrated in time numerically. The momenta are chosen as 
(top) Pt = 0.2 ao, v = 0.4 (JoTo and (bottom) pr ~ 0.06 ag, 
V = 0.2 croTQ. The time axis is scaled as (t/tq)^^^. 



2w„ 



(93) 



dgF{6,6\pT,v) 



g=g'=eo 

-0, (94) 
(95) 



where Cjp±^ — \/p\ + v'^ 1^% + for transverse modes 
and for longitudinal 



modes 



we employ Wp 



Even though the choice of the ini- 
tial conditions is somewhat arbitrary we find that the de- 
pendence on the initial conditions is rather weak as long 
as high momentum modes are sufficiently suppressed to 
avoid cut-off dependencies. In this situation the early- 
time dynamics is dominated by soft modes and the spec- 
tral shape is governed by the primary instability after a 
short period of time. In terms of the original field vari- 
ables and its derivatives with respect to proper time the 



initial conditions (92)-(95) translate to 



(To /6A^ 



3ro V A 



(96) 



for the macroscopic field, where the first derivative origi- 
nates from the time-dependent rescaling of the fields and 
accounts for the decay. Similarly for the fluctuations one 
finds 



F(ro,ro,PT,;^)-^e---/'^o 



(97) 



dr F{t, t',pt,v) 



dr. F{T,T',pT,iy)l=r' 

F{tq,to,pt,v) 
3to 

■F(ro,To,pT,t^) 
9r2 



(98) 



+ ^^-^J^o . (99) 



transverse fluctuations F±^ {t,t,pt,i^)- For times when 
the considered modes are outside the resonance band 
we assumed constant solutions in terms of the confor- 
mal variables. One observes that the evolution is indeed 
well described by the quasi-static approximation. 



APPENDIX D: INITIAL CONDITIONS 



The initial conditions are formulated in terms of the 
conformal fleld variables as 



= 



(92) 



for the macroscopic fleld. We choose the initial fluctu- 
ations to be exponentially suppressed at high momenta 
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